Microscopic Structure of a Vortex Line in a Dilute Superfluid Fermi Gas 
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The microscopic properties of a single vortex in a dilute superfluid Fermi gas at zero temperature 
are examined within the framework of self-consistent Bogoliubov-de Gennes theory. Using only 
physical parameters as input, we study the pair potential, the density, the energy, and the current 
distribution. Comparison of the numerical results with analytical expressions clearly indicates that 
the energy of the vortex is governed by the zero-temperature BCS coherence length. 
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The trapping and cooling of dilute Fermi gases is an 
increasingly active area of research within the field of ul- 
tracold atomic gases. Several experimental groups now 
trap and cool alkali atoms with Fermi statistics reach- 
ing temperatures as low as ~ 0.2Tp with Tp being the 
Fermi temperature of the system [1]. These gases are ap- 
pealing to study due to the large experimental control of 
their properties and the microscopically well-understood 
two-body effective interaction between the atoms. One 
of the main goals of the present experimental effort is to 
observe a phase transition to a superfluid state predicted 
to occur below a certain critical temperature T c [2]. A 
fascinating prospect of such a superfluid is the formation 
of quantized vortices. The combination of theoretical and 
experimental studies of vortices in Bose-Einstein conden- 
sates (BECs) have produced a number of beautiful results 
in recent years [3]. For fermions, the nature of vortices 
in different systems such as type-II superconductors, su- 
perfluid 3 He, and neutron stars is a classic problem with 
a vast literature [4-6] . 

One problem of fundamental importance is to calcu- 
late the energy of a vortex for T = 0. This energy is 
defined as the difference between the energy of the su- 
perfluid state with the vortex present and without. Di- 
viding this value by the angular momentum per parti- 
cle in the vortex state gives the frequency at which the 
vortex becomes the thermodynamic ground state of a ro- 
tating system. It is rather surprising that despite the 
large amount of work concerned with the structure of a 
vortex for a fermionic system, there is no clear result re- 
garding this specific problem. This is in contrast with 
the case of a dilute bosonic superfluid, where the Gross- 
Pitaevskii equation allows an analytical calculation of the 
vortex energy for T = if quantum fluctuations are ne- 
glected [7]. The equivalent theory relevant for fermions, 
Ginzburg-Landau theory, is unfortunately only valid for 
\T — T c \/T c <C 1 making an analytical calculation of the 
energy for T = more complicated. In this letter we 
present the first ab-initio calculation of the vortex en- 
ergy based on microscopic theory. 



We consider in the following a two-component gas of 
neutral fermionic atoms with mass m in a cylinder of ra- 
dius R. We take the number of particles in each spin 
state N a to be the same as this is the optimum situation 
for superfluidity [2]. In the dilute limit, the interaction 
between the atoms in the two spin states a —1,1 can 
be well described by the contact potential g8(r) where 
g = 4Trh 2 a/m and a < is the s-wave scattering length 
describing low energy collisions between | atoms and j 
atoms. In the zero-temperature limit that we are treat- 
ing here there are no intra-component collisions [8]. Re- 
cently, two papers have calculated the T = vortex en- 
ergy for a Fermi superfluid under these assumptions. Us- 
ing phenomenological models, the energy of a unit circu- 
lation vortex was estimated in Ref. [9] to be 
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where n a — kp/Qn 2 is the density in each of the two 
components, and £bcs = Ti 2 kp /nm&o is the BCS co- 
herence length with Ao the bulk value of the superfluid 
gap; BCS theory predicts Ao = 8e _2 /2exp(— 7r/2fc^ \a\) 
[10]. The effective Fermi momentum, kp, is defined as 
H 2 kp/2m = fj, — gn a = ft, where /i is the chemical po- 
tential, such that it includes the effect of the Hartree 
mean-field. It was argued in Ref. [9] that a microscopic 
calculation for T = would yield D to be a constant 
~ 0(1) independent of &f and \a\ since the characteristic 
length-scale of a vortex must be expected to be 0(£bcs)- 
The value of D depends on the phenomenological model 
used: If the vortex is modeled as a cylinder of radius £bcs 
containing a normal stationary fluid, surrounded by a ro- 
tating superfluid one obtains D = 1.36. We refer to this 
simple model as the cylinder model. If Ginzburg-Landau 
theory is applied, we get D = 1.65. 

This conclusion was however disputed in the work of 
Ref. [11]. Here it was argued that the characteristic 
length scale of the vortex is much smaller than £bcs and 
the energy correspondingly higher. This is because the 



structure of the vortex core is determined by the lowest 
lying vortex states. These states are formed out of exci- 
tations around the Fermi level with typical wavelengths 
<~ kp 1 , and following the conclusions based on the analyt- 
ical and numerical solutions of the Bogoliubov-de Gennes 
(BdG) equations [12,13] it was argued that the important 
length scale of the core region is £1 = 4/Trk F \a\ -C £bcs 
in the dilute regime [11]. Using £1 as the size of the 
vortex core, a calculation identical to the one given in 
Ref. [9] leads to a vortex energy given by Eq. (1) but 
with D ~ £bcs/£i 3> 1 in the dilute regime. Thus, the 
energy was predicted to be significantly higher than what 
was estimated in Ref. [9] . Note that D is now not a con- 
stant but depends on kp and \a\. 

It is presently not clear which of the two quite different 
predictions is correct and thus what the energy of the 
vortex actually is. In order to settle this question, we now 
present a microscopic calculation of the vortex energy 
using the assumptions given above. The BdG equations 
describing the superfluid state read [4] 
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with H Q (r) = -h 2 V 2 /2m - \i + U(r), and where U(r) = 
gn a (r) = g^2„ ^(r)! 2 is the Hartree field and A(r) = 
—g u v (r)v*(r) is the gap function. To avoid having to 
introduce an arbitrary high energy cut-off in the theory, 
we have used a zero range pseudo-potential scheme giving 
rise to a regularized coupling constant g, when calculat- 
ing A(r) [14]. The result is a well-defined theory using 
only physical parameters as input. Further details of this 
and the numerical techniques used to solve these equa- 
tions will be given elsewhere. Once the self-consistent 
solution is obtained for a given coupling strength and 
chemical potential, the energy is given by E = (H — pN) , 
where H is the Hamiltonian of the system [4] , and N is 
the number operator. Defining E(r) = 2J2„ \ v v( r )\ 2 
E can be calculated using 



E = -jd 3 r |i?(r) + i[|[/(r)| 2 + |A(r)| 2 ] 



(3) 



We now solve the BdG equations for the gas in a 
cylinder of height L z and radius R 3> £bcs- Since 
£bcs is a decreasing function of kp\a\ we use two dif- 
ferent cylinder sizes. For kp\a\ < 0.43 (fc^H > 0.43) 
we take R = 44.1/im (R = 12.6/im) and L z — 12.5/zm 
(L z — 5/um). We find the lowest energy superfluid state 
of the system by setting A(r) = A(p, z) where p is the 
perpendicular distance from the axis of symmetry of the 
cylinder, and z is the axial coordinate. For a vortex state, 
we assume the form A(r) = exp(— i<p)A(p, z) where <f> 
is the azimuthal angle around the cylinder axis. This 
corresponds to a vortex line with unit circulation along 
the cylinder axis. In both cases U(r) = U(p,z). In the 
vortex-free case, the cylindrical symmetry dictates that 



Cooper pairs form between particles with angular mo- 
mentum hv and — hv along the cylinder axis whereas 
in the vortex case pair constituents have angular mo- 
mentum hv and — h(v +1) [4]. Once the two solutions 
with and without the vortex are obtained, the energy 
per unit length of the vortex line can be determined as 
£ v = (E v — E )/L z where E v denotes the energy of the 
vortex state and Eq the energy of the state without a 
vortex, both obtained from Eq. (3). Throughout the 
succeeding analysis we consider a gas of 6 Li atoms with 
a = —955 ao- The density in each hypcrfinc state is cho- 
sen to be a few times 10 13 cm~ 3 . These parameters are 
appropriate for on-going experiments [15]. The numerical 
solution of the BdG equations is obtained within a Bessel 
function discrete variable representation in p [16,17], and 
periodic boundary conditions along the vortex axis. 
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FIG. 1. The energy density of the superfluid gas. Note 
that on this scale the energies of the system with and without 
a vortex are indistinguishable, as the energy cost associated 
with vortex formation is much smaller than the total energy. 

In Fig. 1, we plot the total energy density E/V of 
the superfluid gas, where V is the volume of the cylin- 
der, and E is given by Eq. (3), as a function of the ef- 
fective interaction strength /cfH- The x's and »'s are 
obtained from a self-consistent numerical solution of the 
BdG equations for two different values of R, whereas the 
line is the ground state energy per unit volume of a bulk 
Fermi superfluid, 
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Here iV(0) = mk F /2n 2 h 2 is the density of states at the 
Fermi level. This expression is obtained by integrating 
analytically Eq. (3) for a homogeneous gas where the 
u(r)'s and w(r)'s are simple plane wave states. The first 
three terms in Eq. (4) give to the energy of a homoge- 
neous gas in the normal phase within the Hartree-Fock 
approximation and the last term is the condensation en- 
ergy. We see that there is good agreement between our 
numerical results and the analytical formula. The slight 
discrepancy is due to boundary effects at the edge of the 



cylinder, where the density of particles vanishes. Note 
that we are at the limit of the weak coupling regime 
kp\a\ -C 1, appropriate for dilute gases. For the pur- 
poses of comparison with analytical results, however, it 
is important to calculate properties for the widest pos- 
sible range of £bcs, subject to the condition £bcs "C R 
which ensures that the gap function can heal to its bulk 
value before becoming suppressed at the cylinder surface. 
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FIG. 2. The energy of the vortex in terms of the param- 
eter D. The dashed and solid lines correspond to the ana- 
lytical predictions of Ref. [9] and [11], respectively, and the 
numerical results are indicated with »'s (R = 12.6pm) and x 's 
(R = 44.1pm) with the average D represented by the dotted 
line. The inset depicts £ v with lines giving the analytical 
prediction of Eq. 1 using D — D. 

In Fig. 2, we plot the numerically calculated energy 
of the vortex £ v for varying kp\a\. To compare with 
the analytical predictions, we parameterize £ v by the 
variable D appearing in Eq. (1). The dashed line cor- 
responds to the prediction D = 1.36 [9] and the solid line 
to D — £bcs/£i [11]. We see that the two predictions for 
D have a completely different dependence on kp\a\. The 
important conclusion is that the numerical results con- 
firm D ~ 0(1) being a constant independent of fcjr|a| in 
agreement with Ref. [9]. On the other hand, the predic- 
tion D = £bcs/£i yields a qualitatively incorrect result. 
We note that the kink in £ v is due to -R/£bcs being differ- 
ent for the two cylinder sizes, whereas the spread in D at 
kp\a\ = 0.43 is indicative of our numerical accuracy. The 
numerical value D ~ 2.5 is higher than the prediction of 
phcnomenological models in Ref. [9] . This is as expected 
since these models only can yield the correct order of 
magnitude of the constant inside the logarithm. Thus, 
the length scale determining the energy of the vortex is 
~ £bcs and not £i . 

To examine this in more detail, we plot in Fig. 3 the 
numerically calculated profile of a vortex for two repre- 
sentative values of fcp|a|. Close to the vortex core, only 
the lowest-energy (bound) states contribute to the or- 
der parameter; these give rise to the observed Friedel 



oscillations, which have a wavelength on the order of 
kp . We see that the length scale defined as £i = 
lim p ^o[A( j o, z)/pA oc ]~ 1 giving the slope of A(r) at the 
vortex core actually is much smaller than £bcs as pre- 
dicted in Ref. [11,12]. Here, Aoo is the value of |A(r)| far 
away from the vortex core with A M ~ Ao as expected. 
However, as the distance p from the vortex core increases, 
the slope decreases and A(r) reaches the value A^ on a 
length scale ~ £bcs and not £i. To quantify this, we 
use the cylinder model of the vortex with a vortex ra- 
dius £2 = x£bcs to calculate £ v . This yields Eq. (1) but 
now with D = (1.36) x jx. The equation D = 2.5 then 
gives x = 0.42. Thus £2 = 0.42£bcs is the length scale 
determining the energy of the vortex. Again, it should 
be emphasized that x ~ 0.42 is a constant over the large 
range of £bcs used in the calculations thereby verifying 
that indeed £bcs determines the length-scale relevant for 
the energy as discussed in [9]. The cylinder model of 
A(p) with the correct radius £2 = 0.42£bcs is plotted in 
Fig. 3. 
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FIG. 3. The vortex profile A(p)/Aoo for two values of 
&;fH. The thick solid line corresponds to 25,000 atoms per 
hyperfine state and a transition temperature of 0.045/iK while 
the thick dashed curve is for N a = 66, 500 giving T c = 0.23pK. 
For both curves R = 12.6pm. The thin solid and dashed 
lines depict the cylinder model of the vortex with radius 
£2 = 0.42£bcs for the two fci?|o| values. The inset shows 
the full &fH = 0.43 solution. 

To examine the supcrfluid flow associated with the vor- 
tex giving rise to the angular momentum, we plot the 
current density j s (p) given by 
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in Fig. 4 for fcp|a| = 0.59. Because the normal component 
carries no current by construction, the total current den- 
sity may be written as j s (p) = 2n s v s with the supcrfluid 
velocity w s — e r j,h/2mp and the supcrfluid (or atom-pair) 
density n s . We plot n s (p) defined in this way in Fig. 4. 
Note that unlike dilute interacting Bose gases at zero 



temperature, the supcrfluid density in a dilute Fermi su- 
pcrfluid does not have the same behavior as the order pa- 
rameter A(r). As expected, n s (p) ~ n a far away from the 
vortex core and n s (p) — > for p — > 0. The dotted line in 
Fig. 4 gives the analytical result j<p(p) — hkp/(6w 2 mp), 
which agrees well with the numerics away from the vor- 
tex core. As a self-consistency check on the numerics, 
the angular momentum per particle along z is found to 
be exactly h/2, corresponding to one unit of angular mo- 
mentum ft, per Cooper pair as expected. 
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FIG. 4. Upper panel: the current density j s (p) = jV(p) e ^ 
(solid line) and its asymptotic form (dotted line). Lower 
panel: the density n a (p), the superfluid density n s (p) and 
the gap function A(p) normalized to their theoretical bulk 
values. In both panels fcj?|o| = 0.59. 

As shown in Fig. 4, the presence of a vortex in a Fermi 
system does not lead to any significant change in the par- 
ticle density n a (p) [4], in contrast with a dilute Bose gas, 
where the density is minuscule on the vortex line [18]. 
Direct observation of the vortex core (now commonplace 
for BECs) is therefore not likely. The quantized currents, 
and therefore the presence of superfluidity, can be readily 
detected using at least three approaches, however. One 
of these is the collective mode spectrum. When no vortex 
is present, excitations carrying equal and opposite angu- 
lar momentum along the z-axis are degenerate in energy. 
The vortex currents lift this degeneracy since the rota- 
tional symmetry is removed. The resultant splitting of 
the surface modes is proportional to the angular momen- 
tum of the gas [9,19]. This technique has been used to 
infer the presence of a vortex in a trapped BEC [20]. A 
second approach was demonstrated in a recent experi- 
ment where the precession rate of the scissors oscillation 
mode was used to measure the quantized angular mo- 
mentum per particle with great accuracy [21]. A third 
method is spatially selective Bragg scattering; the super- 
fluid currents modify the Bragg momentum conservation 
conditions, giving rise to a strongly anisotropic outcou- 
pled atomic beam [22]. 

In conclusion, we have accurately determined the en- 
ergy of an isolated vortex line in a dilute superfluid Fermi 



gas at zero temperature. The results clearly indicate that 
the BCS coherence length sets the scale for the vortex en- 
ergy and therefore the critical frequency for vortex sta- 
bility. 

We acknowledge valuable discussions with 0. Elgar0y, 
C. J. Pethick, and B. I. Schneider. 
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